Inflationary models inducing non-Gaussian metric fluctuations 
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(Dated: February 2, 2008) 

We construct explicit models of multi-field inflation in which the primordial metric fluctuations 
do not necessarily obey Gaussian statistics. These models are realizations of mechanisms in which 
non-Gaussianity is first generated by a light scalar field and then transferred into curvature fluctu- 
ations. The probability distribution functions of the metric perturbation at the end of inflation are 
computed. This provides a guideline for designing strategies to search for non-Gaussian signals in 
future CMB and large scale structure surveys. 



Inflation generically predicts Gaussian initial metric 
fluctuations with an almost scale invariant power spec- 
trum |Q. It is interesting to determine how general this 
property is. With the advent of large scale structure 
and CMB surveys it will be possible to test in details 
whether the initial conditions were Gaussian or not and 
in the latter case how they can be characterized. At 
the moment, no non-Gaussian signal has been detected 
in CMB data [|| but the number of modes that can be 
probed is still small. In large-scale structure surveys 
the number of independent modes that are observed is 
large but the difficulty arises from the non-linear gravita- 
tional dynamics fl which can shadow the primordial non- 
Gaussianity [|] . Having at our disposal models of infla- 
tion in which non-Gaussian adiabatic metric fluctuations 
are generated can then serve as a guideline for designing 
strategies for detecting primordial non-Gaussianities. 

In single field inflation, the slow-roll conditions, if valid 
throughout the period during which the seeds of the large 
scale structures are generated, prevent the generation of 
observable primordial non-Gaussianities. The reason is 
that the potential needs to be both flat enough for the 
fluctuations to develop and steep enough for the non- 
linearity to be significant. Note however that starting 
from a non-vacuum initial state || or allowing for sharp 
features in the potential could lead to primordial non- 
Gaussianities. In those cases these features will be local- 
ized in a narrow band of wavelengths and will affect the 
shape of the density fluctuation power spectrum. 

The situation is richer when more than one light scalar 
field are present during inflation. In multi-field inflation 
models, a mixture of adiabatic and isocurvature fluctua- 
tions [Q , that can be correlated or not |j] , are generically 
produced. It opens the door to the generation of primor- 
dial non-Gaussianity Q because the non-linear couplings 
can be much stronger in the isocurvature direction than 
in the (inflaton) adiabatic direction [l(|. When the 
fields are coupled, exchange between isocurvature and 



adiabatic modes can be important Jll| which makes it 
possible to generate non-Gaussianity on a large band of 
wavelengths while preserving an adiabatic slow-roll type 
primordial power spectrum ]12] . 

In a previous study |p^| , we have presented the build- 
ing blocks for such models. It is based on the genera- 
tion of non-Gaussian isocurvature fluctuations which are 
subsequently transferred to the adiabatic mode due to a 
bend in the classical inflaton trajectory. The necessary 
ingredients for such models are (i) the existence of a light 
scalar field that has non-linear coupling, typically quar- 
tic, and (ii) a coupling term in the potential leading to 
isocurvature-adiabatic mode mixing. This latter mech- 
anism can be given a simple geometrical interpretation. 
On super Hubble scales, curvature fluctuations are in- 
duced by fluctuations of the total expansion in different 
parts of the universe. The relative duration of their in- 
flationary phase can be affected if inflaton trajectory is 
bent. Fluctuations in the isocurvature modes generate a 
bundle of parallel trajectories along which inflation lasts 
longer or shorter depending on whether they lie in the 
outer or inner part of the bent (see Fig. 1 of Ref. Q). 

In Ref. we characterized the expected statistical 
properties of the metric fluctuations that were shown to 
be the superposition of a Gaussian and a non-Gaussian 
contributions of the same variance. The relative weight 
of the two contributions is related to the bending angle 
in field space. We explicitly computed the probability 
distribution function (PDF) of the non-Gaussian contri- 
bution which appears to be described by a single new 
parameter, [see Eq. (||) below]. These results were ob- 
tained without the use of any explicit model. This letter 
is dedicated to the building of this kind of inflationary 
models involving either two or three fields. The shape of 
their PDF will be computed and compared to our previ- 
ous analytical expression . 

In general the curvature fluctuations, 1Z, can be com- 
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puted on superhorizon scales as |Tj 
K = HSt. 



(1) 



Isocurvature fluctuations can then induce curvature fluc- 
tuations only if in some way or another they can change 
the relative total number of e-fold. 

This can actually be obtained in a very simple way in 
hybrid type model of inflation with three fields. In that 
case, one field, <fi, is the inflaton; the second field is a light 
scalar, Xi with quartic coupling A < 1 and the third field, 
a, is coupled to the two others so that the end of inflation 
is triggered when a undergoes a phase transition. To be 
more explicit, let us examine the following model 



V((/>,x,cr) = 
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--a {<p cos a + x sin ex.) 
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where cto is the final vev of tr and a is the mixing angle 
between <j> and x in their coupling to a. 

Inflation ends when the effective mass of a vanishes, 
that is when 



g ((f> cos a + xsina) — 2^(Tq = 0. 
The value of 6 at the end of inflation is 



cos a 



(3) 
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For tfi > (f>endi cr = and the two fields evolve inde- 
pendently: <f> drives the inflation while x develops non- 
Gaussianity. The amount of non-Gaussianity of x then 
depends only on A and on the total number of e-fold be- 
tween horizon crossing and the end of inflation. 

When a is non-zero, fluctuations of x induce curva- 
ture fluctuations because they change the time at which 
the phase transition occurs. Thus the x-hiduced curva- 
ture fluctuations read (assuming H is basically constant 
during the inflationary period) , 



K ~ H5t e 



3H 2 sin a 



cos a 



(5) 



We present on Fig. [I] the result of a numerical inves- 
tigation in which the super-Hubble evolution of the x 
modes are explicitly computed. The Klein-Gordon equa- 
tions for the two fields and the Friedmann equations are 
integrated for a set of initial conditions leading to a bun- 
dle of classical trajectories. The parameters of this ex- 
ample were taken to be A = 1CP 2 and an evolution for 
a number of e-folds equal to N e = 50. The shape of the 
PDF is compared to our analytical expression |12[| , 



P (x)d X = \ 77- 




exp 



(6) 




FIG. 1: The shape of the probability distribution function of 
the isocurvature modes in case of a quartic coupling, ^| (solid 
line) compared to the analytic PDF of Eq. (^|) (dashed line) 
and to Gaussian distribution (dotted line). In this example, 



A = 10~ z and N e 



50. 



which depends only on the variance of x an d on the pa- 
rameter z/3 quantifying the amount of non-Gaussianity. 
The variance is directly proportional to the value of the 
Hubble constant at horizon crossing, H c . The parameter 
P3 is proportional to A and is explicitly given by 



1/3 = -XN e /{'6H 2 c ) w -0.60. 



(7) 



On Fig. 0, the numerical results are compared to the 
expression (|^) and to a Gaussian distribution of the same 
variance. It confirms that the analytic form (||) is a very 
good approximation of the PDF. The value of ^3 in this 
plot is jjjl vz — — 0.15 /a 2 ,, which implies that the value 
of the kurtosis of distribution is 



s 4 



(x 4 )/(x 2 ) 2 -3 = 4^ 



(8) 



The curvature fluctuation PDF is related to the PDF of 
the field %. In the limit where the inflation ends abruptly 
the relation is straightforward, Eq. (||). It is more com- 
plex in general and depends on the details of the phase 
of inflation. 

The parameters of the potential (|^) are tuned, as usual, 
so that the amplitude of the fluctuations are compati- 
ble with the observations. No further fine-tunings are 
required for the initial conditions: whatever the initial 
value for the field %, it rolls down rapidly towards x = 
to follow the valley bottom, as described in the previous 
paragraphs. 

The model (^) reads as a simple polynomial potential 
of order four but involves three fields. It is actually possi- 
ble to build models involving only two fields. In this case, 
it is clear that the x trajectory cannot be straight and 
has to bend. A possible explicit model can be obtained 
when the zero-mode trajectory is given by 



Xo(t) = Xoo tanh 



acj> (t) 



Xc 



(9) 



where a and Xoo are free parameters. A model to get such 
a trajectory can be constructed by considering a term of 
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the form (x — Xo) 4 m the potential, as first advocated in 
Rcf. p2[, so that the potential might take the form 



V{4>,x) = \m 2 2 + ^ 
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(10) 



Although it gives the adequate trajectory (|9|) in field 
space, it is non-polynomial and does not induce signif- 
icant metric non-Gaussianities. The reason is that the 
X fluctuations are squashed as soon as the trajectory is 
bent, the bundle of trajectories being drastically focused 
much before any significant bending. This is a generic 
features of all models in which the bending of the trajec- 
tory is driven by the non-linear coupling term. 

A way to construct a model in which the x fluctuations 
are not squashed is to consider the potential 



V(ct>, X ) 
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Now, the zero mode trajectory (</>o, Xo) is not determined 
by the non-linear coupling. At lowest order in the cou- 
pling parameter A and in the slow-roll regime, the Klein- 
Gordon equations for <j) and x can be integrated and lead 
to the trajectory (j^). 

In Fig. ^, we depict the behavior of the super- horizon cj) 
and x fluctuations (left panels) and the curvature fluctu- 
ations that are induced by the initial x fluctuations (right 
panels). The parameters of the model for the simulation 
of Fig. U were chosen to be 



m = 10~ 7 - Mpi. 
Xoo = 1.5 Mpi. 



a = 1 

A = 5 x 10~ 5 , 



(12) 
(13) 



so that it gives realistic fluctuation amplitudes and e- 
fold number. The initial conditions are set up at 
</>init = 4Mpi.. At this time the coupling constant is 
then about 0.15 and the number of e- folds during the 
inflationary period is about 72. We integrate the Klein- 
Gordon equations and the Friedmann equation for the 
set of classical trajectories starting from different ini- 
tial conditions around (^>i n ;t , 0) Gaussian distributed with 

((^init) 2 )-^ 2 )-^. 



As can be observed on the left panels the joint evolu- 
tion of the </> and X fluctuations is quite complex. This is 
due to the fact that the fluctuation distribution is shaped 
by the nonlinear couplings during the trajectory bend. 
The resulting curvature fluctuations can however be sim- 
ply described. They are given by the sum of those in- 
duced by the <p fluctuations, 7^.^, (Gaussian distributed) 
and those induced by the x fluctuations, 1Z X . Both distri- 
butions have approximately the same variance. For the 
parameters (|l2|- 13|) , our numerical computations show 
that actually the r.m.s. of the former is about 18iJ c , 
the r.m.s. of the latter 15 H c . The right panels show 
the PDF shape of the non-Gaussian component of the 
curvature fluctuations, 1Z X . They follow the distribution 



function given in Eq. (|6|) as long as | S4 1 < 0.7 at least for 
less than 3<r events. The complete PDF of the curvature 
fluctuations can then be obtained as the convolution of 
the PDFs of the two components, IZ^, and 1Z X . 



In this letter, we have introduced a series of multi- 
field inflationary models leading to the generation of non- 
Gaussian primordial curvature fluctuations. We have 
been able to build explicit models with either two or three 
scalar fields. In the case of two-field models, it is difficult 
to design a simple potential: it has to be constructed 
in such a way that the non-linear term does not damp 
the curvature fluctuations. In the case of three fields, 
we have shown that models as simple as polynomial po- 
tentials of order four work. In both cases we have illus- 
trated our models by numerical computations and have 
shown that the analytic result (||) derived in Ref. |l^] is 
a very good approximation to the primordial curvature 
fluctuation PDF when non-Gaussianities in the x dis- 
tribution remain modest. Let us emphasize that in the 
two-field models the potential has to be tuned but that 
this is not the case for the three-field model. Moreover, 
the initial conditions do not require to be fine-tuned, as 
explained before. This description of the primordial non- 
Gaussianity is the starting point of further investigations 
on its various observational aspects, both for CMB and 
large scale structure surveys fuf. 
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Note that the variance of the field value evolves due to 
the non-zero effective mass of %■ 



